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We report on the first next-to-leading BFKL study of the cross section and azimuthal decorrel- 
lation of Mueller Navelet jets. This includes next-to-leading corrections to the Green's function 
as well as next-to-leading corrections to the Mueller Navelet vertices. The obtained results for 
standard observables proposed for studies of Mueller Navelet jets show that both sources of cor- 
rections are of equal and big importance for final magnitude and final behavior of observables, 
in particular for the LHC kinematics investigated here in detail. The astonishing conclusion of 
our analysis is that the observables obtained within the complete next-lo-leading order BFKL 
framework of the present contribution are quite similar to the same observables obtained within 
next-to-leading logarithm DGLAP type treatment. The only noticeable difference is the ratio the 
azimuthal angular moments (cos 2^) / (cos (p) which still differs in both treatments. 
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1. Introduction 

The high energy regime of QCD is one of the key questions of particle physics. In the semi- 
hard regime of a scattering process in which s » —t, logarithms of the type [ailn(5/|?|)]" have 
to be resummed, giving the leading logarithmic (LL) BaUtsky-Fadin-Kuraev-Lipatov (BFKL) [1] 
Pomeron contribution to the gluon Green's function describing the ^-channel exchange. To reveal 
this effect, various tests have been proposed in inclusive [2], semi-inclusive [3] and exclusive pro- 
cesses [4]. The basic idea is to select specific observables minimizing usual coUinear logarithmic 
effects a la DGLAP [5] with respect to the BFKL one: the involved transverse scales should thus 
be of similar order of magnitude. We here consider the Mueller Navelet jets [6] in hadron-hadron 
colUders, defined as being separated by a large relative rapidity, while having two similar transverse 
energies. In a DGLAP scenario, an almost back-to-back emission is expected, while the allowed 
BFKL emission of partons between these two jets leads in principle to a larger cross-section, with 
a reduced azimuthal correlation between them. We review results of Ref. [7] where both the NLL 
Green function [8] and the NLL result for the jet vertices [9] are taken into account. 



2. NLL calculation 

The two hadrons collide at a center of mass energy producing two very forward jets, whose 
transverse momenta are labeled by Euclidean two dimensional vectors ky i and ky a, while their 
azimuthal angles are noted as <j>j^i and <j>j^2- The jet rapidities yy i and are related to the longi- 
tudinal momentum fractions of the jets viaxy = We restrict ouselves to fixed rapidities and 
transverse momenta. For large xj^i and xj^2, colUnear factorization leads to 

i,/"^, ^ = E/'dxi C dx2Ux,)Ux2) ... , (2.1) 

d|k7,i|d|k7,2|dy/,idy7,2 ^io io d|k7,i|d|k7,2|dy/,idy7,2 

where /a,b are the parton distribution functions (PDFs) of a parton a (b) in the according proton. 
The resummation of logarithmically enhanced contributions are included through A:r-factorization: 

|.|.'''^|' = / d(^'y,i#/,2 / d2kid2k2Va(-ki,xi)G(ki,k2,^)Vb(k2,;c2), (2.2) 

d|ky,i|d|k7,2|dy/,idy7,2 J J 

where the BFKL Green's function G depends ons = xiX2S. The jet vertices ^ were calculated at 
NLL order in Ref. [9]. Combining the PDFs with the jet vertices one writes 

— = / ^J,l^J,2 /d2kid2k2C&(ky,i,Xy,i,-ki)G(ki,k2,5)<I>(kj,2,^7,2,k2), 

d|k7,i|d|ky,2|dy7,idy7,2 J J 

where ^(k7,2,X7,2,k2) = y^dx2/(x2)V(k2,X2). (2.3) 

In view of the azimuthal decorrelation we want to investigate, we define the coefficients 

'^m=j #7,2 COS (m((/>y,i - (^7,2 ^^^i d^k2 <I>(ky,i , -ki )G(ki , k2, j)<I>(ky,2,^y,2 , k2) , 

from which one can easily obtain the differential cross section and azimuthal decorrelation as 
d|k,,i|d|kt'|dy.,idy,,2 =^" '"'^ (cos(m(p)) = (cos(m(,;..,i-,^>.,2-^))) = ^. (2.4) 
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The guiding principle of the calculation is then to use the LL-BFKL eigenfunctions 

£„,v(ki) = -L(k2)'^-3^-^,^ (2.5) 
7rv2 

although they strictly speaking do not diagonalize the NLL BFKL kernel. In the LL approximation, 
^^ = (4-35^,o)y dvC,„,v(|k/,i|,^y,i)Qv(|k7,2|,x/,2)f -j , (2.6) 

where Cm.v{\kj\,xj) = J d^jd^kdxf{x)V {k,x)E,„^v(k)cos{m^j) , (2.7) 
and (0{n,v) = Nca,/7iXo(\n\,^ + iv) , mth Xo{n,Y) = 2^{l) (y+ 'i) (l - Y+ ^) . The 



master formulae of the LL calculation ( p^q , [2.7[ ) will also be used for the NLL calculation, the 
eigenvalue now turning to an operator containing a V derivative, which acts on the impact factors 
and effectively leads to a contribution to the eigenvalue which depends on the impact factors. 

At NLL, the jet vertices are intimately dependent on the jet algorithm [9]. We here use the 
cone algorithm. At NLL, one should also pay attention to the choice of scale ^o- We find the choice 

2 

of scale So = .JsoJsoJ with 5o,i = j rather natural, since it does not depend on the momenta 
ki 2 to be integrated out. Besides, the dependence with respect to of the whole amplitude can 
be studied, when taking account the fact that both the NLL BFKL Green function and the vertex 
functions are dependent. In order to study the effect of possible coUinear improvement [10], 
we have, in a separate study, implemented for « = the scheme 3 of the first paper of Ref. [10]. 
This is only required by the Green function since we could show by a numerical study that the jet 
vertices are free of 7 poles and thus do not call for any collinear improvement. In practice, the use 



of Eqs. ( 2.6 , 2.7 ) leads to the possibility to calculate for a limited number of m the coefficients C„. 



as universal grids in v, instead of using a two-dimensional grid in k space. We use MSTW 2008 
PDFs [11] and a two-loop strong coupling with a scale /Ir = ylK/jPlky^ . In order to compare our 
analysis with DGLAP NLO approaches [12] obtained through the NLL-DGLAP partonic generator 
Duet [13], for which symmetric configurations lead to instabilities, we here display our results for 
|k/_i I = 35GeV, |k72| = 50GeV (see Ref. [7] for symmetric configurations). 



3. Results 

Fig. |l|a and [T]b respectively display the cross-section and the azimuthal correlation as a func- 
tion of the relative jet rapidity Y, for the LHC design center of mass energy ^/s = 14TeV. This 
explicitely shows the dramatic effect of the NLL vertex corrections, of the same order as the one 
for the Green function [14]. In particular, the decorrelation based on our full NLL analysis is 
very small, similar to the one based on NLO DGLAP. The main source of uncertainties is due to 
the renormalization scale Hr and to the energy scale ^/sq. This is particularly important for the 
azimuthal correlation, which, when including a collinear improved Green's function, may exceed 
1 for small }JLr = pLp. The only remaining observable for which a noticeable difference can be 
expected between BFKL and DGLAP type of treatment is the ratio (cos 29) / (cos <p) . 

The NLL analysis presented here could be extended to describe forward jets production. 

Work supported in part by the French-Polish Scientific Collaboration Agreement POLONIUM, 
the grant ANR-06-JCJC-0084 and by a PRIN grant (MIUR, Italy). 
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Figure 1: Differential cross section (a), azimuthal correlation (cos(p) (b) and ratio (cos2(p)/(cos(fi) (c) in 
dependence on Y for |kj 1 1 = 35GeV, |ky.2| = 50GeV. The errors due to the Monte Carlo integration are 
given as error bands. Blue: pure LL result; Brown: pure NLL result; Green: combination of LL vertices 
with the collinear improved NLL Green's function; Red: full NLL vertices with the collinear improved NLL 
Green's function. Dots show the results of Ref. [12] obtained with DiJET [13]. 
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